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Abstract 



We study the novel decays B* ^ 7 7 and D* ^ D 7 7 using the frame- 
work of the Heavy Meson Chiral Lagrangian (HMxL) to leading order in chiral 
perturbation theory. The branching ratios of these decays are expressed in 
terms of the strong gB*(D*)B{D)-K aiid the electromagnetic gB*{D*)B{D)'y cou- 
plings, thus providing a possible tool for their determination. In the charm 
case, using the experimentally determined ratios DTr)/(D*^'~^ 
1)7), we are able to express the branching ratio as a function of the strong 
coupling only. We thus find 1.6 x 10"^ < Bt{D*° ^ 77 ) < 3.3 x 10"^ 
for 0.25 < g < 1, where g is the strong coupling of HM^L. In the beauty 
sector, the Bt{B*^ — > 5*^77) which we estimate to be in the 10^^ — 10~^ 
range is a function of both gs* Bn and gB*B'y Its behaviour does not afford an 
unambiguous determination of these couplings except for the region of high 
g values like g > 0.6. The expected two-photon differential distributions are 
presented for both B*^ — > i?'^ 7 7 and D*^ D^j'j, for different values of 
the couplings involved. 
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I. INTRODUCTION 



The heavy vector mesons B* and D* (of spin-parity 1~) decay via spin-flip electromag- 
netic or strong interactions to the well studied pseudoscalar ground states B and D. The 
decays of D* are known to proceed either as a strong transition D* Dtt with a final pion 
with momentum of about 40 MeV or as an electromagnetic one — > D 7 , with a final 
photon with momentum of about 140MeV. The situation is different for B* which has a 
mass of 5324.9 ± 1.8 MeV; since the mass difference Mb* — Mb is only 45.8 MeV, there is no 
strong B* decay and the radiative process — > B^is the dominant decay mode for B*. 

In the present paper we study another possible electromagnetic decay, the two-photon 
decay processes B* —>■ 7 7 and D* —>■ D 7 7 which were not considered previously in the 
literature. In addition to the intrinsic interest in these novel modes, we point out that their 
study could provide information on the strong couplings gs^Bw and Qd^d-k and on the elec- 
tromagnetic ones gB'B-y, go'D-y- The strong couplings are directly related to the basic strong 
coupling g of the effective heavy meson chiral Lagrangian, which describes the interactions 
of heavy mesons with low-momentum pions. 

There is a major dissimilarity between the possibility of measuring the couplings in the 
charm and beauty sectors, as a result of the different mass difference between the respec- 
tive vector and pseudoscalar mesons. In the charm sector, the experimentally measured 
branching ratios of the D*^, D*^ decays into the allowed Dt: and D7 modes lead to re- 
lations between gD*D'K and gD*D-y Henceforth, the D* — > D77decay under study here is 
expressible in terms of the strong coupling only and would provide a convenient tool for its 
measurement. 

On the other hand, in the beauty sector where the B* ^ Bn decay is forbidden by phase 
space, the gB*B-K coupling is not directly accessible. And although the B* S 7 decay is 
experimentally detected, the direct measurement of its strength is an unlikely proposition 
at present, in view of the smallness of the expected value of its decay width. 

However, as we show in this paper, the two photon decay B* — > i? 7 7 branching ra- 
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tio turns out to be a function of gB^Bw and Qb^b-^, which opens the possibihty for their 
determination, especially in particularly favorable regions of the [qb'B-ki Qb'B'^] parameter 
space. 

Our analysis sing les out the neutral modes B*^{D*^) B^{D^) 7 7 as the more rele- 
vant ones in relation with the determination of the the couplings under consideration. In 
addition to the direct radiative transition which is the sole contribution in neutral decays, 
and on which we concentrate in this paper, there is also the two-photon decay arising from 
bremsstrahlung in the charged B*^[D*^) B^[D^)'^ channel. Since this radiation over- 
whelms the direct mode, as we will show, one has to resort to the investigation of the neutral 
modes if one aims for a cleaner determination of the strong couplings. 

In section II we review the present status of the main decays of B* and D* ^ with which 
the rare two photon decays must be compared. In section III we present the theoretical 
framework of our approach. Section IV contains the explicit treatment of the decay ampli- 
tudes. In the last section we summarize our predictions and we discuss certain features of 
the calculation. 

II. THE EXPERIMENTAL AND THEORETICAL STATUS OF B* AND D* 

PRINCIPAL DECAYS 

The vector mesons B* were firstly observed by the CUSB collaboration at the Cor- 
nell Electron Storage Ring (CESR) by detecting the photon signal from the radiative decay 
B* ^ B J . This signal of 46 MeV photons was confirmed in improved CUSB-II measure- 
ments, with the vector mesons produced at CESR at the T{5S) resonance 0. Recently the 
process B* i?7has been observed also at LEP with the various detectors in a sample 
of over 4 milion hadronic decays. The rate of the i?*-meson production relative to that 
of B mesons is found generally to be consistent with the expectation from spin counting. It 
is expected that this production rate will be maintained in future B experiments at hadron 
machines, like BTeV at Fermilab and LHC at CERN, where samples of the order of 10^*^ -B's 
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are expected. Then, a fairly high sensitivity can be achieved in the study of B* decays and 
measurements of B* branching ratios of the order 10^^ — 10^*^ could be accessible. 

In view of the small mass difference AM{B* - B) = 45.78 ± 0.35 MeV g|, which forbids 
strong B* decays,the electromagnetic transition B* —>■ i? 7 appears as the main decay of B*. 
This decay has been studied in a variety of theoretical models, including quark models , 
the chiral bag model P| followed by effective chiral Lagrangian approaches for heavy and light 
mesons P HTll] , potential models [|T2| , [T3| and QCD sum rules 0,1^. The predictions of these 
calculations span a range of nearly one order of magnitude for the expected decay widths, 
between r(5*°(5*+) ^ 5°(fi+)7) = 0.04(0.10) KeV m and T{B*\B*+) B°{B+)-f) = 
0.28(0.62) KeV |^, with most of the calculations HJ7|,p|, p!5|] giving values closer to the larger 
values of Ref . IB! . 



The D* meson was discovered more than twenty years ago [0 and has been subsequently 
studied in several experiments at different accelerators e.g. |T7|-p2|. The D*^ (M = 2010.0 ± 
0.5 MeV) and D*^ (M = 2006. 7±0. 5 MeV) have relatively little phase space for strong decay 
into D + 71. The current PDG averages |^ for the measured branching ratios of the observed 
decays are Bt{D*+ D+ti^) : Br(D*+ D^n+) : Br(D*+ D+7) = (30.6 ± 2.5)% : 
(68.3 ± 1.4)% : (1.1 + 2.1 - 0.7)% and Br(L)*° ^ DV) : Br(D*° ^ D°-f) = (61.9 ± 2.9)% : 
(38.1±2.9)%. The most recent experiment on D*'^ decays |^T| gives more accurate branching 
ratios for the three decay channels as follow (30.7 ± 0.7)% : (67.6 ± 0.9)% : (1.7 ± 0.6)%. 

Although there are, by now, good data on the branching ratios, there is still no absolute 
measurement of any of the partial decay widths. The tightest upper limit has been estab- 



hshed by the ACCMOR Collaboration at CERN [||] from the measurement of 127 
events using a high-resolution silicon vertex detector, to be r{D*^) < 131 KeV. The other 
closest limit, obtained by the HRS collaboration |]18|, gives upper limits of 1.1 MeV and 2.1 
MeV for the total decay widths of the charged and neutral D*' s. 

The decays of D* have also been treated extensively in a plethora of theoretical models. 
Many of the papers we mentioned concerning the B* decay p|-p!^] discuss also the D* decays. 
In addition, we want to mention the early approaches p3ip4| with effective Lagrangian 



including symmetry-breaking, a relativistic quark model |2^, the study of D* decays using 
the chiral-bag model which contains pion exchange effects (pion loops) ||26[, the use of QCD 



sum rules p^, a chiral model with Mc — > oo |28| and the comprehensive analysis of Kamal 



and Xu ||2^. Recently [Q, the strength of the various decay channels of D* has been 
extracted from an analysis of the experimental branching ratios by the use of the chiral 
perturbation theory. 

In the D* case, the theoretical calculations again span an order of magnitude range for 
the prediction of the absolute decay widths, from a small width of T(D*^) ~ (3 — 10) KeV 
|1|,|28| to r(D*°) ~ (60 - 120) KeV ij2|j2|], including fairly large uncertainties. It should 



be emphasized at this point that the chiral bag calculation (x) has offered the best estimation 
for the branching ratios |2§ Br^(D*+ D+ti^) : Bi^{D*+ D^n+) : Br^(D*+ D+-f) = 
31.2% : 67.5% : 1.3% and Br'^(D*° ^ D^n^) : Br^(D*° ^ D^) = 64.3% : 35.7%. 

The recent experiments have confirmed these relative ratios and dispersed the 

puzzling features which prevailed previously concerning the radiative branching ratio and 
the relative ratios of D*^ strong channels (see, e.g. Ref. [^). The prediction of the chiral 
bag model P§] is T{D*+ all) = 79 KeV, T{D*^ all) = 59.4 KeV. Several of the other 



calculations result in fairly similar values [p|,|5|,pU| as well as predicting for T(D*^) a value 
approximately 25% smaller than for T{D*^). There are also calculations in which these 



widths are nearly equal |lTT| , p^ or, on the contrary, calculations giving T(D* ) to be at least 
twice larger than T{D*°) p|,|7|,|8 . 



The experimental and theoretical survey we presented here is obviously of direct relevance 
to our calculation as the absolute values of B* and D* widths will affect the probability of 
the B* — > 57 7 and D* 1)7 7 future detection. 

III. A MODEL FOR TWO-PHOTON TRANSITION. 

Substantial progress has been made in recent years in the treatment of the interactions 
of heavy mesons containing a single heavy quark with low momentum pions, by the use of an 
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effective Lagrangian |]31|-p3|, the so-called "Heavy Meson Chiral Lagrangian" (HMx L), which 
embodies two principal symmetries of Quantum Chromodynamics (for a comprehensive 
review of this theoretical framework and its applications, see [§4|)- leading order 

in an 1/Mh expansion {Mh is the mass of the heavy meson) and the chiral limit for the 
light quarks {mi 0,1 = u,d,s), the Lagrangian carries flavor and spin symmetry in the 
heavy meson sector, as well as SU{3)l ® SU{3)ii chiral invariance in the light meson one. 
We adopt this framework for the calculation of the processes we study here, namely the 
B*{D*) B{D)'y'y decays, and we shall use it to display the possible usefulness of these 
transitions for the determination of the Qb'B-k , gD*D-K and gs'S-y couplings. 

The heavy vector {B* or D*) and pseudoscalar {B or D) mesons are represented by a 
4x4 Dirac matrix H , with one spinor index for the heavy quark and the second one for the 
light degree of freedom. 

1 + 



H 



B*r-B^, , H = ^oHho, (3.1) 



where v is the meson velocity, v'^B* = and B*, B are the respective annihilation 
operators of the meson fields. We shall usually refer henceforth to B* i? 7 7 with the 
understanding that the same treatment holds for D* ^ D 7 7 . However we shall specify 
the two channels separately when numerical or other specific features make it necessary. 

The relevant interaction term of the HM^L, representing the coupling of heavy mesons 
to an odd number of pions, is given by ||3l| , |32[| 

^i?MxL = {Ha7,75A^abH,) (3.2) 
where the axial current is 

A^='- {eo^^ - ^d^^^) (3.3) 

and ^ = exp(zA^ / /) , with A4 being the usual 3x3 matrix describing the octet of pseudoscalar 
Nambu-Goldstone bosons. The axial coupling constant g is one of the basic parameters of 
HM^L, which is of direct import to our problem, a, b denote light quark flavours (a, 6=1,2,3) 
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and f is the pion decay constant, / = 132 MeV. Expanding the axial current and using the 
first term = —{l/f)d^Ai + .... one obtains the effective Lagrangian representing B*B- 
pion and i?*i?*-pion interactions, which are the relevant ones in our problem. Thus, 



''-eff 



■'^B*d^'MB^ + h.c. 
f ^ 



(3.4) 



The dimensionless B*B'k coupling is defined as p4 | 



(3.5) 



where is the polarization vector of -B*, with the physical coupling given by the limit 
m^. We use the same normalization convention as in [p4| . 



Throughout this work, we assume that the variation of Qe'E-k with in the region of 
our treatment can be safely neglected. Likewise, we define 



(7r(g)5*(^;i,ei)|5*(i;2,e2)) = gB*B*.{q')eap^.Ae^2q''v'i 



(3.6) 



with the same remarks as above. 

We also note that the isospin symmetry requires 



93*317 = g3* + 30T. 



-93'-03+7T- 



-V2g_ 



3* + B+TvO 



^93*°BOnO 



(3.7) 



and the gB*B-K so defined is the commonly used in the literature. 

The (|377| ) relation holds similarly for gB*B*n, gD*DTT and gD*D*-K couplings. 
Now, from ( p.2|) -( pl^ , using the definition ( p.5|) one has 



93*371 



gD*Dn 



2Mb 
2Md 



g — g3*B*TT 



f 



-g — gD*D*TT- 



(3.8) 



Note that, in deriving (|3.8|) one assumes B, B*{D, D*) mass degeneracy. In order to calculate 



the B* ^ i? 7 7 decay width we use the interaction Lagrangian (|3.4| ) to the leading order 
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in chiral perturbation theory, which is an appropriate tool here, in view of the smallness of 
Mb* - Mb. 

The calculation of the radiative processes B*{D*) B{D)'~f'~f obviously requires the 
incorporation of the electromagnetic interaction in our Lagrangian ( p.4| ), which is per- 
formed pD| , |55| by the usual procedure of gauging the Lagrangian with the U{1) photon 
field. This leads to the replacement of the derivative operators in the Lagrangian by co- 



variant derivatives containing the photon field, explicitly exhibited in [^. Nevertheless, 
the new Lagrangian still does not provide for couplings to induce the observed B* ^ -B 7 , 
D* — > D 7 magnetic dipole transitions. This necessitates the introduction of an additional 
term in the Lagrangian, a contact gauge invariant interaction proportional to the electro- 
magnetic field Ffj^i,, which is given by [p0| , p5| 



^ = ^Tr (^Ha(J^yF^yHi,5ab 



(3.9) 



where is the strength of this anomalous magnetic dipole interaction, having mass dimension 
[1/M]. 

Additional terms arising from an 1/Mh expansion exist PD|J53[ ], however several of them, 
including the radiation of the heavy quark can be absorbed in equation ( ^.91 ) by redefining 
In the present paper, we shall consider ^ as an effective coupling, representing the strength 
of the B*{D*)B{D)'y transition. 

Expanding H in terms of the components, one obtains for the additional electromagnetic 
interaction 



Cls = -e/iF^'^ tB;^B: + t^p^,v''{B^Bl + h.c.) 



a/ Rt R* 



(3.10) 



which exhibits B*B'~f and B*B*'~f couplings with equal strength, as given by the heavy 
quark symmetry. From ( p. 101 ) we obtain the respective vertices, which are 



{-f{k, e)B*{vi, ei)\B* {v2, €2)) = e/iMs. (ei.A;e.e2 - e2.ke.e1] 
{^{k,e)B{v,)\B*{v2,e2)) = -leMB*^ie,,^pe^k''v^ei 



(3.11) 
(3.12) 
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The propagator of the heavy vector meson is given by —i{g'^'^ — v^^v")/2 [{v.k) — A/4] , 
where A = Mb* — Mb and f , k are the velocity and the residual momentum. The propagator 
of the heavy pseudoscalar meson is z/2 [(f ./c) + 3 A/4] ||34|| . 

Now, considering Lagrangians ( |3.4| ) and ( |3.10| ), as well as the axial anomaly responsible 
for the 7r° — 77, we classify the diagrams contributing io B* ^ 577 in the leading order 
of chiral perturbation theory as follows: 

There is the diagram B*^ i?'^"7r" 5^77 (Fig.l), via a virtual pion (the some- 
what different situation in D* D77will be analyzed in the last section). This diagram 
contains the known strength of the pion axial anomaly. Then, there is the loop graph 
B*^ (i?*^7r^) — s> B^'-yy (Fig.2), with the photons radiated from the virtual charged pion 
in the loop, with additional graphs of the same class, as specified in the next section. The first 
graph is proportional to the gs^Bn coupling, while the loop graph contains the gB*B-KgB*B*ir 
product. In addition we have the tree level diagram with two insertions of the magnetic 
operator defined in ( |3.10| ), leading to B*^ B*^y B^jj (Fig.3). This graph depends 
only on the /i magnetic moment. Finally we have a class of one loop diagrams which involve 
both the magnetic moment and the strong coupling, which is exhibited in Fig. (4-6). We did 
not include the contribution of diagrams containing three heavy meson propagators which 
is negligible. Needless to say, the determination of gB*B-K would be simpler, should the first 
two graphs dominate. However, this is not true for D* decay, while it can be true for the 
B* decay for an opportune range of parameters as we discuss in the next section. 

We remark at this point that corrections to ( p.2| ) which arise from higher terms in the 
1/Mh expansion as well as in chiral breaking have also been investigated p0|j3^j3^j37[] . A 
comprehensive inclusion of these corrections in the calculation of the two-photon decays of 
heavy vector mesons is beyond our scope in this first treatment of these processes. Neverthe- 
less, we note that we shall use physical masses for the degenerate doublet of heavy mesons 
both in the loop propagators and in the decay calculations, moreover the chiral loops in- 
cluded are themselves of order 1/Mh- The terms we include are the leading ones in chiral 
perturbation theory, and are of the same order in an 1/Nc expansion; moreover to this order 
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there are no counterterms [0S 



It is appropriate to mention now that the gB^B-K, gD*Dn couphngs were estimated in 



recent years by the use of a variety of theoretical techniques, hke QCD sum rules |T^,p7 



soft pion approximation and other methods P, p!0| , |30| , |36| , |37| , ^ . Generally, the values 
of g obtained in these works are in the range g = 0.25 — 0.7, significantly smaller than 



the quark model result of = 1 or of modified quark models P^j41|j4^ which brought this 
value slightly below one. The most recent determinations of g include an analysis |Q of 
various theoretical approaches which leads to a "best estimate" of g = 0.38, a recent lattice 



determination giving 0.42(4) (8) and the analysis of Stewart |Q which incorporates 
symmetry breaking terms in the Lagrangian and obtains g = 0.27ig 4. 

Finally, the experimental limit 131 KeV p2[ puts an upper limit oi g < 0.71, 

using r{D*+ D^Ti+ + D+TfO) = {g"^ / A-k p)\p^\^ . 

The existing theoretical estimates we mentioned, give 0.04KeV < T{B*^ ^ B^-f) < 
IKeV and O.lOKeV < T{B*^ l) < IKeV, where we allowed for a slightly higher 

upper limit. We redefine the magnetic coupling in eq.( p.lOD to a dimensionless quantity 
/i = M^./i = gB*OBO^, and /i+ = MB*fi+ = gB*+B+-yy then the above limits give the following 
ranges, 2.2 < < 11.0, and 3.5 < fL+ < 11.0. 



IV. THE DECAY AMPLITUDES. 

We present now the explicit expressions of the decay amplitudes, which in our approach, 
to leading order in chiral perturbation theory, consist of the contribution of the anomaly 
graph (Fig. 1), the tree level graph (Fig. 3) and the loop graphs (Figs. 2,4-6). 

In presenting the differential decay distribution, we use the following variables: 

t={p-hf = {p' + hY 

u={p-k2f = {p' + hf (4.1) 
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with 



t + s + u = ML + M\ 



(4.2) 



where fci, ^2 are the four momenta of the two photons and p,p' are the four-momenta of the 
decaying B* and the final B respectively. The allowed ranges for s and t are 



< s < {Mb" 



Msf, t- <t<t^ 



{Ml, +Ml-s)± J {Ml, + Ml- sf - 4Ml,Ml 



(4.3) 



The amplitudes are given for B*'^ ^ 5*^ 7 7 and we shall remark on the changes appearing 
in D° 7 7 whenever required. The amplitude from the anomaly graph, mediated by 

a pion is: 



A 



anomaly 



{B* 



V271J 



(4.4) 



where eB*,Gi,€2 are the polarization vectors of the heavy vector meson B* and the two 
photons respectively. There are additional contributions from rj, rj' which are not specified 
in ( |4.4| ). As we shall describe in the next section, their contribution is rather small and we 
may safely neglect them at this stage. 
For the tree level graph we find 

47ra/i^ 



At,ee(5*"^ 50 77: 



M: 



B 



t-Ml, 



u-Ml, 



(4.5) 



The loop contribution which depends only on the strong coupling is given by a sum of 
several diagrams. In addition to that explicitly shown in Fig. 2 there are diagrams with one 
photon radiated by the virtual pion in the loop and the other emitted from the B*B*7t, 
B*B'K vertices or both photons emitted from these vertices, or both photons emitted from 
the loop by the 7r7r77 vertex. In the limit of photon momenta small compared to the pion 
mass, which we find to be a suitable approximation, the class of diagrams of Fig. 2 gives 
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OigB*BTTgB*_B*-K 

SttMr* 



X 



(4.6) 



Finally other loop contributions to the B* 7 7 decay come from diagrams where 
both the strong coupling and the magnetic one are involved. There are diagrams with one 
photon radiated by the virtual pion in the loop and the other emitted from the ingoing 
B* particle through the B*B*[B)'y vertices (Fig. 4), or from the outgoing particle B* which 
becomes B through the B*B^ vertex (Fig. 5), or from the B in the loop which becomes B* 
through the B*Bj vertex (Fig. 6). 

The amplitude corresponding to Fig.4 with a B*B*j vertex is: 



g lla{MB* - MB)gB*B*ngB*Bnl^ 

^77) = 7^3772 eaa7<5 



167rM|* 



+ {k^{eB*.e2)-e'i{eB*.k2)) 



t-Ml, 
{p - k2yelk1 
u- Ml, 



(4.7) 



The amplitude corresponding to Fig.4 where a B*B^ vertex replaces the B*B*^ appear- 
ing in Fig.4, is 



t - Ml, 



ei^k2^ej 



u 



Ml, 



The amplitude corresponding to Fig.5 with B in the loop is: 



(4.8) 



Afoop(5*°- 5°77) 



2,amig%,B^H 



87rMl, 



{p - k^rei^kl^ej , {p-k2re^,k'(el 



t - Ml, 



+ 



u - Ml, 



and the one with B* in the loop is 

*° ^ 5O77) 



Afoop(^ 



\QttM^ ^0/37(5 ^B*'^ ^-qaip^ ixvariyg V-j- -r g^V ) 



(p-k,)^(p-k,yele^k^^ , (p-k2)Hp-k2re1elk'(' 



t-Ml, 



+ 



u - Ml, 



(4.9) 



(4.10) 
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Finally the amplitude corresponding to Fig. 6 is: 



(4.11) 



In the B* decay, the pions are the sole contributions in the loop, while in the 
D* ^ D'-f'-f calculation we include both pions and kaons. 
Let us call A the sum of all the amplitudes: 

where Aioops is the sum of all the amplitudes ^loops) = 1, 6 which come from the loops. 

The square of the above amplitude, when averaged over the initial spin and summed 
over the final spins is: 

\A\' = It.\A\'^ (4-13) 

spins 

where we have included the factor | in order to take into account two identical particles 
in the final state. 

The differential decay rate of photon energy is obtained by integrating the following 
expression over the variable t, 

\A\^dt. (4.14) 



ds (27r)3 32M|, 

There is a major difference in the anomaly contribution of the B* and D* decays. Since 
the 7r° appears in the physical region in the D*^ 7r° decay we have to isolate the on- 

shell tt" decay in the D*^ D° 7 7 mode. Hence, for the D* case we limit ourselves in the 
integration of dV/ds to a region which goes from s = up to 20 MeV away from the pion 
mass. 

Using now the physical masses of Mb* , Mb, Md^Md |@] and the eq.( p.8|) we obtain for 
the decay rates of 5* — -B 7 7 , the following expression: 
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T{B* 577) = 3.40 X IQ-'V + 1-53 x lO'^V + 4.81 x IQ-'V/i 

+ 1.53 X lO^^V/i + 4.71 X 10'^^ gjf + 9.81 x lO^^V/i^ 
+ 2.65 X lO^^^^V + 1-38 X IQ-^^Jx^ + 2.67 x 10~^V/W+ 
+ 2.90 X 10"^VAi+ + 8.94 x 10-^^g^fijx+ 

+ 9.11 X 10"20^VAi+ + 7.21 X 10"^V/i+- (4.15) 

As we can see this is a function of both the strong couphng g and the magnetic dipole strength 
which represent the effective 5'b*oso^ and gB*+B+-y couphngs. For the D* 1^7 7 we 

have: 

V{D* D 7 7 ) = 2.52 X 10"^ V + 5.85 x lO^^V + 1-79 x 10"^ V/i 

+ 4.43 X IQ-^^g^Ji + 1.30 x lO^^^^^u^ + 3.50 x 10""^V 
+ 2.42 X 10~^^^V + 2.18 X 10-^V + 1-01 x 10-"^^/i+ 
+ 2.95 X 10~^VAi+ + 2.11 X 10-^VAi/^+ 

+ 1.70 X 10"^^^VAi+ + 2.05 X 10"^VAi+- (4-16) 

In this case we can relate the magnetic couphng to the strong couphng using the existing 
experimental informations on V{D*^ D^n^) : T{D*^ D^-f) of (61.9 ± 2.9)% : (38.1 ± 
2.9)%, which gives rise to the relation: fi ~ 6.6^, and T{D*+ ^ D° 7r+ ) : r(D*+ -> D+7) 
of (67.6 ± 0.9)% : (1.7 ± 0.6)%, which gives rise to the relation: /i+ ~ 1-7(7. Then we can 
write the decay width solely as a function of g, which is a crucial step in the engagement of 
this decay as a tool for measuring g. 

T{D* D-f-f) = 2.52 X 10"^ V + 5.66 x 10"^ V 

+ 4.76 X lO^y + 3.64 x 10"^°^^ + 1.53 x lO^V- (4.17) 



We used in eqs.( |4.l5| )-( |4.16| ) the same notation of fi, /i+ for the magnetic moment strength 



in both the charmed and beauty sectors although they are probably not equal for the physical 
processes. However since in the charm case, the magnetic coupling has been related to the 
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strong one, the /i, /i+ will denote in the rest of the paper the strength of the B^'~^'y 
transitions. 

The experimentally measured branching ratios of D* — D n , D* ^ D'j lead to relations 
between /i, and g modulo an unknown phase. We have allowed also for the possibilities 
of negative relative signs among /2, /i+ and g and as it turns out, this affects only slightly the 
numerical picture, due to the fact that the main contribution is given by quadratic terms. 
In eg. ( [4. IS] ) we give for comparison the expression obtained for fl = —6.6(7,/!+ = —1.7g. 



r{D* I^7 7)^=neg = 2.52 X 10""/ + 5.66 x lO""/ 

+ 4.70 X lO^V - 3.64 X lO'^^g^ + 1.53 x 10"V- (4.18) 

This ambiguity will be further discussed in the next section. The difference between the 
rates of B* and D* is mainly due to the different phase space. 

In discussing the two photon radiative decays, we shall refer in the next section to the 
following quantities 



BRiB* ^ 577) = .7.1 = .o^. V (4.19) 



BR(D ^ D77) = ^^^;o^ = r(Z^.o^^o^) + r(Z^.o_^o,ov (4-20) 



V. DISCUSSION AND CONCLUSIONS. 

The formalism we have presented refers to the decays of the neutral heavy vector 
mesons B*^, D*^, as it will the numerical analysis of our results, which will be given be- 
low. For the charged decays, B*^ —>■ i?+7 7and D*^ D^'j'-f one has to consider also 
the bremsstrahlung emission which appears in diagrams of Fig. 3 and Fig.4 and additional 
ones. The bremsstrahlung radiation comes from the initial or the final charged particles. To 
give an idea of this effect we have calculated the part of the bremsstrahlung amplitude which 
is due to radiation from the final B~^ particle in the amplitude B*^ {B^l) -8^77. It is 
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^brem = ^TTagB^B^e^^e^e'lp'' 



{t-MD ^ {u-MD 



(5.1) 



This amplitude (and the ones given by B*^ radiation) have to be added in order to get the 
full amplitude for the B*^ B"^ 7 7 decay. An estimate of the bremsstrahlung decay width, 
from (|5.1| ) only, using for the unknown gB*B'y vertex a value leading to r{B*^ l) = 

0.14 KeV g leads to a decay width of ~ 10"^ GeV for ki,k2 > 10 MeV, considerably 
larger than ([4.15|) . The use of the charged B*^ i?+7 7thus involves a different type of 



analysis in view of the relative size of the different components of \A{B*^ 5+77)^ and 
is less useful for a determination of g. A similar situation is encountered for D*^ 7 7 . 

Thus, we concentrate here on the "safer" neutral decays and we relegate the discussion of the 
charged decays to a separate work, in which we consider the usefulness of T{B*^ B^ 77 ) 



for the determination of gB*+B+-t Ell 



We proceed now to analyse the results on the two decays separately and we start with 
£)*o _^ 7 7 transition for which the rate ( [4.17|) was obtained. 



The many theoretical estimates for g we mentioned in Sec. 3 are spread over the range 
0.25 < g <1 (we also remind the reader that the experimental result ||2^ on the upper limit 
of V{D*^ all) can be interpreted as g < 0.71). Using this range and, eqs. (|4.171 ),( [4.18|) we 
can establish the expectation 

r(D*° ^ DO 7 7 ) ~ (0.022 - 6.73)eV. (5.2) 

The most promising feature of the present analysis arises when we use the existing experi- 
mental informations on V{D*^ ti^ ) : T{D*° D^j) of (61.9 ±2.9)% : (38.1 ±2.9)% 
to transform eq. ([4.20|) into a ratio of T{D*^ 7 7 ) to the total D*'^ width which becomes 



proportional to g^. 
Using 



12n P 



r(D*o ^ /}0^0) ^ ^|p;|3 = 1,25 X 10-VGeV (5.3) 



and the (61.9 ± 2.9)% : (38.1 ± 2.9)% relative branching ratio, we arrive at 
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T{D*° all) = (2.02 ± 0.12) x lO'^^^GeV. (5.4) 



Thus from ( [5.4| ) with ( [1.17| ) we can obtain a branching ratio which depends on g"^ only: 



(0.025 + 0.057^ + 4.76^2 ^ g.SG^^ ^ 1.53^4^ ^ ^g^V 
2.02 X 10-V 



Br(D*° ^ 7 7 ) = • ono, .n-/. ' ^- (5-5) 



With our model for D*^ 7 7 , the measurement of this ratio will thus constitute a 

measurement of the g coupling. Using again the accepted expectation of 0.25 < (7 < 1, we 
predict 

Br(D* ^ D77) = ^ rp,) = (0.16 - 3.3) x 10"^ (5.6) 

A few remarks are in order. Firstly, the sign question The observed branching ratios 
do not afford to estabilish experimentally the sign of g/g£)*oijo^. On the other hand, there 
is theoretical support from the analysis of Stewart on the positive sign of this ratio. 
However even if we assume opposite sign for various pairs of the couplings, we found that 
the changes are rather small, and this is explicitly exhibited in Table 2, and included in 
Eq.(ig). 

The differential distribution in the s-variable can also be used to learn about the value 
of g, due to the fact that the different contributions depend on different powers of g. Finally 
we remark that the contributions from the diagrams exhibited in Fig.4 are rather small, as 
a result of two heavy propagators. The main contributions are those of the anomaly and of 
the graphs of Fig. (2) and (3). In Figs. (7) and (8) we present the differential distributions 
in s for g = 0.7 and g = 0.25. In the latter, the contributions containing a higher power 
of g are diminished and the effect of the anomaly becomes visible in the higher end of the 
spectrum. 

Turning now to the B* ^ B^'-f , we have a rather different situation. Firstly, there 
is only one major decay of B*, namely B* ^ B'-f , which precludes an analysis like in D* 
decays. The branching ratio ( [4.19| ) depends on three parameters, gB*Bn (or g), (7^.0^0^ (or 
/i) and gB*+B+f (o^ /^^)- At this point, we rely on the theoretical estimates presented in 
Section 2 and contain our analysis to the regions given by existing calculations. 
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Now, an inspection of eq. (f4.16|) shows that which appears only in diagram of Fig. 6 



has a very httle effect on the rate, whether /i, /i+ are at the lowest or at the highest end of 
their value, for any value of g. Hence, we continue our analysis in the parameter space of 
[g,^] only. 

In table 2 we present the values of Br(i?* -B 7 7 ) for different values of g and the two 
extreme values of /x, corresponding to T{B* ^ B 'y) = 40eV and IKeV. Again, assuming 
that relative negative signs are possible we give in the last column the Br for /i = —2.2. 
Clearly, a branching ratio in the 10^'' — 10^^ range will not allow to pinpoint accurate values 
for the two couplings. 

Nevertheless if the branching ratio turns out to be in the 10~^ range, it can only be 
caused by large values of g, say g > 0.6. 

We wish also to mention an additional scenario: the g coupling will probably be measured 
directly in D* decays, or indirectly from Br(D*° D'^'yy) or other methods. With this 
knowledge, Br(i?*° B^'j'j) becomes a function of gs'^BO-y only and it could provide the 
desirable measurement of this coupling. This is a very interesting issue, since as pointed 
out already some time ago |]2B|, there is no other possibility of measuring the width of the 



5* — i> i? 7 decay with presently known techniques. 

In Figs. (9), (10) and (11) we give the differential distribution of dr{B* B'y'y)/ds for 
g = 0.5 and three different values of fi. Clearly, once g is known one may use accurate 
differential distributions to distinguish between different fi values. 

At this point we also wish to make some remarks on the similar decays of the strange 
heavy vector mesons, B*^ — * -8^77 and Ds*~^ Df'j'j which were not mentioned so 
far. In both these cases the pion anomaly is further suppressed, since 

7r° can proceed only by isospin violation, e.g. via rj — n^ mixing. On the other hand, both 
decays can proceed via chiral loops with charged i^-mesons in the loops, 5*° ^ {K+B*-) 
5° 7 7 and D*j (K^D*^) -D+77. However, one must add that for D*+ ^ £)+77 
there is the complication of the bremstrahlung and we shall disregard it here. 

We calculated therefore only the B*^ decay, and the situation is quite similar 



to that encountered in B* decay; therefore we do not repeat this analysis here. 

Before concluding we comment on a few points which were neglected in our treatment. 

1 We calculated also the contribution to the anomaly term of a virtual 1] exchange for 
the D* —>■ 1^7 7 decay. The inclusion of rj modifies our result in eq.( [4.16[ ) by a factor 
of 1 + ,^^^ . Since go'Dr) and go'D-w are expected to be comparable, this is a small 
effect. 

2 We neglected the off-shell g^-dependence of the anomaly which could have some effect, 
especially in the D*-decay. This should be included in a more detailed treatment. 

To summarize, we have used the Heavy Meson Chiral Lagrangian to present the first 
treatment of the rare B*^ 5*^77, D*^ 7 7 decays. The decay rates depend on 

the strong gB*B-K 1 Qd'Dtt couplings and on the strength of the magnetic dipole transitions 
gB*B-/, gD*D'y The strong couplings are expressed in the chiral Lagrangian by the strong axial 
coupling g. 

We have shown that Br(D*° 7 7 ) can be given as a function of g only, and as such, 

it would provide an appropriate tool for its measurement. For the conventionally envisaged 
range 0.25 < ^ < 1 we calculated 1.6 x 10"^ < Bt{D*° D°77) < 3.3 x 10^^ 

On the other hand, Br(i?*° 5^77) is a function of gB*BTT , gB'^B^-y 9B*+B+'y The 
latter coupling has little effect on the branching ratio. Nevertheless, one cannot determine 
specific values for the first two couplings from the measured branching ratio, unless g is in 
the higher part of its expected range, say 0.6 — 1. 

The differential dV/ds distributions in both cases can be used as additional help for 
extracting the values of the coupling constants. If the value of gs^Bn will turn out to 
be in the "measurable" range, it will be of great interest to check the HMxL relation 

gB*B7T/gD*Dn = Mb/Md. 
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TABLE I. Predictions for the various D* — > D'yy decay for various values of g (eqs.(5.9) and 
(5.12)) 



9 BR(£>*)(/Z = 6.6ff,/Z+ = 1.7ff) BR{D*){fi = -6.6g,fi+ = -1.7g) 

g = 0.25 1.7x10"^ 

g = 0.38 3.9 x 10^^ 

g = 0.5 6.9 X 10-6 

g = 0.7 1.4 X 10-^ 

g = l 3.3 X 10-5 



TABLE II. Predictions for the various B* B'yy decay different values of g and n 



9 


BR{B*){j2 = 2.2) 


BR(S*)(/i = 11.0) 


BR(S*)(/i = -2.2) 


g = 0.25 


3.1 X 10-7 


1.9 X 10-6 


2.4 X 10-'^ 


g = 0.38 


1.3 X 10-6 


2.2 X 10-6 


9.4 X lO-'^ 


5 = 0.5 


3.7 X 10-6 


2.9 X 10-6 


2.7 X 10-6 


3 = 0.7 


1.4 X 10-5 


5.5 X 10-6 


9.2 X 10-6 


3 = 1 


4.8 X 10-5 


9.0 X 10-6 


3.8 X 10-5 




k2 



B*(D*) B(D) 
FIG. 1. The anomaly graph for B*{D*) B{D)jj. 
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1.7 X 10-6 
3.7 X 10-6 
6.3 X 10-6 
1.3 X 10-5 
2.9 X 10-5 




+ 



B'(D') B'(D') B(D) 

FIG. 2. Loop graph for B*(D*) B(D)'~f'y which depends only on the strong couphng. Addi- 
tional graphs of this kind are discussed in the text. 



k, 



B*(D*) B*(D*) B(D) 



FIG. 3. Tree graph for B*{D*) B{D)'y'y. 



+ crossed graph 



7r(7r,K) ~l~ crossed graph 



B*(D*) B*(D*) B*(D*) B(D) 

FIG. 4. Loop graph for B*{D*) —>■ B{D)'j'y wich involve both the magnetic and the strong 
couplings. Instead of B*{D*) in the first propagator we can have B{D). 



fci 



fe 



7r(7r, K) 



+ crossed graph 



B*(D*) B*(D*) B*(D*) B(D) 

FIG. 5. Loop graph for B*{D*) B{D)j^ wich involve both the magnetic and the strong 
couplings. Instead of B*{D*) in the loop we can have B{D). 
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fcl 



7r(7r, K) 



+ crossed graph 



B*(D*)B(D) /B*(D*) B(D) 



FIG. 6. Loop graph for B*{D*) — > B{D)^^ wich involve both the magnetic and the strong 
couplings. 
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2 4 6 8 10 12 

s GeV'-2) 

FIG. 7. The differential decay width dT{D* — > D'^^)/ds (eV) as a function of s with the value 
g = 0.7. 
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s (10"-3 GeV"2) 

FIG. 8. The differential decay width dT{D* — > D'y^)/ds (eV) as a function of s with the value 



g = 0.25. 
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0.5 1 1.5 2 

s (lO-'-a GeV"2) 

FIG. 10. The differential decay widths dT{B* B'yj)/ds (eV) as a function of s, with the 
value g = 0.5 and /i = 5.7. 
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